Experimental and theoretical analysis of the upper critical field in FSF trilayers 
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The upper critical magnetic field Hc2 in thin-film FSF trilayers is studied experimentally and 
theoretically in geometries perpendicular and parallel to the heterostructure surface. The series 
of samples with variable thicknesses dpi of the bottom and dF2 of the top Cu4iNi59 F-layers are 
prepared in a single run, utilizing a wedge deposition technique. The critical field Hc2 is measured in 
the temperature range 0.4 — 8 K and for magnetic fields up to 9 Tesla. In order to properly interpret 
the experimental data, we develop a quasiclassical theory, enabling one to evaluate the temperature 
dependence of the critical field for an arbitrary set of the system parameters. A good agreement 
between our experimental data and theoretical predictions is demonstrated for all samples, using a 
single set of fit parameters. 

PACS numbers: 74.25. Dw, 74.45. -l-c, 75.75.Cd 



I. INTRODUCTION 

The upper critical magnetic field Hc2 of an isotropic 
type-II superconductor generally obeys a linear temper- 
ature dependence in the vicinity of the superconducting 
transition temperature Tg^. Deviations from the linear 
T-dependence of Hc2 (T) are usually ascribed to inhomo- 
geneities distributed in the sample, which can broaden 
the resistive transitions R{T) and R{H}^-i. The temper- 
ature dependence of Hc2 is also known to be sensitive 
to the orientation of the magnetic field, if the supercon- 
ducting nucleus size becomes comparable to the char- 
acteristic dimensions of the structure (see, for example, 
Rcfs. 013). In particular, artificially prepared metallic 
multilayers (ML) consisting of alternating superconduct- 
ing (S) and normal metal (N), or of S and insulating (I) 
layers, or even of two different superconductors, S and 
S', show nonlinear Hc2{T) dependences^. 

Among a variety of layered superconducting systems, 
superconductor-ferromagnet (S/F) metallic hybrids at- 
tract special attention because of rich physics of these 
objects^i"— as well as promising perspectives for their ap- 
plications in superconducting spintronicsiSrJ^. Based on 
the unique properties of S/F hybrids, several kinds of 
device physics were proposed such as proximity-effect su- 
perconducting spin-switchingi^"— and Josephson current 
switchingi2ri2£ actuated by an external magnetic field. 
Very recently, magnetic-field-controlled superconductiv- 
ity switching in the S/F proximity systems^lr— and the 
Josephson current switching in the S/F/S junctions^S 
were demonstrated experimentally. 

The influence of the magnetic field on superconductiv- 
ity in S/F hybrid structures was pointed out at the very 
beginning of their studies?^-. It was demonstrated that in 
order to evaluate correctly their critical magnetic fields 
it is necessary to account for the magnetic field penetra- 
tion into both, the superconducting and ferromagnetic 



layers. An approximate single mode anal ysis for the crit- 
ical field^^ was further developed in Ref. [2^ where rigor- 
ous solution of the quasiclassical Usadel equations^^ was 
worked out both for SF-bilayers and for SF-multilayers. 
These studies of the critical fields are important not only 
because of possible influence of the pivoting field on the 
onset of superconductivity in the system, but, in view 
of the large number of physical parameters (from 6 to 
12, see below), also to reduce ambiguities in fitting pro- 
cedures by including more experimental data from the 
same set of samples. 

Being interested in superconducting spin-valve physics, 
proposed by the Beasley group!-- and developed later 
on in Refs. fl^ - flTl we present here a detailed anal- 
ysis of the temperature dependence for the critical 
magnetic field Hc2 of FSF superconducting spin-valve 
cores. We proceed both experimentally, by studying 
Cu4iNi59/Nb/Cu4iNi59 FSF trilayers, and theoretically, 
by evaluating the critical fields with the aid of the qua- 
siclassical formalism of Usadel equations. 

Our paper is organized as follows. In Section II we 
present details of our fabrication procedure of FSF tri- 
layer samples and describe our measurements of their 
temperature dependent critical magnetic field Hc2{T). 
Section III is devoted to theoretical analysis of the prob- 
lem within the framework of quasiclassical Usadel equa- 
tions. Section IV presents a detailed comparison between 
our experimental data and theoretical calculations as well 
as a brief discussion of our key observations. 



II. EXPERIMENT 

A. Thin film deposition and sample preparation 

The three-layer Cu4iNi5g/Nb/Cu4iNi59 structures 
were grown on commercial silicon substrate utilizing a 
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magnetron sputtering machine by the Leybold Company, 
model Z-400. Three mounted targets allowed fabrication 
of the whole multi-structured samples within one cycle 
of sputtering without breaking vacuum in the chamber. 
The interfaces between the constituent thin films were 
found extremely clean^*?. Since our study requires a set of 
specimens with identical thickness of the superconduct- 
ing layer of Nb and varied thickness of the ferromagnetic 
layers, the method of growing wedge-shaped films was 
applied. It allowed us to obtain a set of up to 40 speci- 
mens prepared in a single run on a long substrate under 
identical vacuum conditions. 

Actually, a commercial silicon wafer was cut into pieces 
of size 7 X 80 mm^, cleaned and put into the vacuum 
chamber. After that the chamber was evacuated with a 
turbo-molecular pumping system to a base pressure of 
3 X 10^^ mbar, and the sputtering was done in the at- 
mosphere of argon (99.999% purity) with the pressure of 
8 X 10^^ mbar. At the beginning, each target (75 mm 
in diameter) was presputtered during about 10 minutes 
in order to clean its surface from possible oxides and ab- 
sorbed gases. 

The first Cu4iNi5g layer was grown on an amorphous 
silicon buffer film deposited just before, to isolate the 
structure from gases absorbed by a natural oxide on 
a surface of the substrate. To provide superior homo- 
geneity of the Si and Nb film thickness over the size 
of the substrate the "spray" deposition technique was 
employed'^^''^^. To realize this technique a magnetron 
was moved along the substrate by a custom motorized 
setup providing a uniform coverage of the substrate sur- 
face with the sputtered material. 

To deposit a wedge-shaped Cu4iNi5g layer, the tar- 
get was positioned with the symmetry axis just above 
the substrate edge. A wedge shape of the ferromagnetic 
alloy was obtained due to the native gradient of sput- 
tering aside of the symmetry axis. The deposition rate 
on the thick end of the wedge was about 3-4 nm/sec. 
As in the earlier experiments22i2^, we operated the mag- 
netron in the ac regime to keep the concentration of nickel 
about 59% in the deposited CuNi alloy films. As we al- 
ready pointed out, the superconducting niobium layer 
was grown with the Nb target moving along the sub- 
strate. In this way the reduced effective deposition rate 
of about 1.3 nm/sec was achieved^SiM. Subsequently, the 
second wedge-shaped ferromagnetic layer was sputtered 
on top of the Nb film with the thickness of the both 
Cu4iNi5g layers varied from 2 to 50 nm (more accurate 
data are presented in Sec. II). In order to protect our 
FSF structure from oxidation in ambient conditions it 
was capped by a thin silicon layer. The sketch of the 
whole FSF structure is presented in Fig. [TJ 

To obtain a series of FSF strips with varying Cu4iNi59 
layer thicknesses dp to be used in our Hc2{T,dF) mea- 
surements, samples of equal width (about 2.5 mm) were 
sequentally cut perpendicularly to the wedge gradient 
and correspondingly numbered resulting in the FSF3 
batch. Then, aluminum wires of 50 fim in diameter were 




FIG. 1: (Color online) The sketch of our trilayer FSF sample 
with top (CuNi-T), bottom (CuNi-B) Cu4iNi59 ferromagnet 
layers and superconducting (S) Nb layer. 

attached to the strips by ultrasonic bonding for four- 
probe resistance measurements. 

B. Thickness and composition characterization 

The thickness of the each layer in our FSF structure 
was determined by means of Rutherford Backscattering 
Spectrometry (RBS). He++ ions were accelerated to en- 
ergies of 3.5 MeV by a tandem accelerator. In order 
to avoid channeling effects the samples were tilted az- 
imuthally by 7°, and the backscattered ions were detected 
at the angle of 170° with respect to the incident beam. 

The resulting RBS spectrum is displayed in Fig. [2] 
The thickness of the Nb layer was found to be in the 
ten nanometers range. Further processing was necessary 
to determine the F-bottom and F-top layer thicknesses 
because of the overlap between the Cu and Ni signals 
from the bottom and top CuNi alloy layers. As a first 
step, the Si cap layer signal was identified as the peak 
around the energy of ^-^1980 keV in the RBS-spectrum. 
Due to the energy loss of alpha-particles, backscattered 
from the top CuNi alloy layer and passed through the Si- 
film, the right-hand-side overlapped-peak feature of the 
integral CuNi alloy signal can be identified as Cu sig- 
nal from the top Cu4iNi59 layer. The alpha-particles, 
backscattered from the bottom CuNi alloy layer, loose 
more energy passing through the Nb, top CuNi and cap 
Si-films. That is why, the left-hand-side feature of the 
combined CuNi bottom-|-top signal can be identified as 
a Ni signal from the bottom CuNi layer. The intense 
peak at about 2640 keV results from the combined sig- 
nal of the top Ni and bottom Cu constituents overlapped 
at the same range of energies. Unfortunately, it was not 
possible to determine accurately the concentration x in 
the Cui_a:Nia; alloy, however, our previous experience 
demonstrates^^i^i^i^ that rf-sputtering of the CuNi- 
alloy usually provides x = 0.59 (±1 %). This concen- 
tration is also assumed here within our RBS evaluation. 

The results of the RBS evaluation are presented in 
Fig. [3l Depending on the sample, the graph shows lin- 
early increasing thicknesses of both copper-nickel layers 
from about 1 to 48 nm. The niobium layer thickness is 
nearly constant of around 15.5 nm (dashed line in Fig.|3]) 
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FIG. 2: (Color online) The RBS-spectrum of a single sample. 
The specimen has the following composition: Si-substrate/ Si- 
buffer/CuNi/Nb/CuNi/Si-Cap. 



FIG. 4: (Color online) Resistive transitions for the sample 
#26 of the batch FSF3 in the presence of the magnetic field. 
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FIG. 3: (Color online) Thicknesses of the layers of FSF3 batch 
evaluated by RBS spectrum over the whole strip-like speci- 
men. 



with slight tendency to increase for largest sample num- 
bers. The thicknesses of both CuNi alloy layers were 
found to be nearly the same only for intermediate range 
of the samples (approximately for Nr. 20-^27). How- 
ever, the deviation of the CuNi layers thicknesses over 
the mean value mainly did not exceed 10 %. 



C. Critical fields measurements 

The critical field measurements were performed in a 
^He cryostat equipped with a superconducting solenoid 
providing magnetic fields up to 17 T. The temperature 
was controlled in the range of 0.4 — 8 K with an accuracy 
of about 1 mK. Resistivity measurements were performed 
in Oxford Instruments '^He-'^He Heliox insert, using the 
conventional four-probe method. In order to avoid ther- 
moelectric voltages, alternating of the current polarity 



was performed during our measurements. The resistive 
transitions for the sample #26 are displayed in Fig. |3]for 
different values of the magnetic field. The critical tem- 
peratures Tc assigned to the applied field, i.e. the upper 
critical magnetic field Hc2 , were determined as midpoints 
of the resistive transitions R{T). The critical values of 
the magnetic field were measured for two field orienta- 
tions - perpendicular and parallel to the sample plane, 
H^(T) and i/j,'2(T), respectively. Both the critical tem- 
perature values and the slope of the curves (T) in the 
vicinity of Tc are listed in Table ID 

Figure [S] displays the temperature dependencies for 
perpendicular (left panel) and parallel (right panel) crit- 
ical magnetic fields for the FSF3 series with the niobium 
layer thickness value c^Nb = 15.5 nm. One observes that 
both critical fields demonstrate a non-monotonous de- 
pendence on the ferromagnetic Cu4iNi5g layer thickness. 
This observation requires an explanation which will be 
obtained by means of our theoretical analysis developed 
in the next section. 



III. THEORY 

In this Section we elaborate a theoretical analysis of 
the superconducting phase transition in F1SF2 trilayers 
in the presence of an external magnetic field. A sketch of 
our trilayer is shown in Fig. [S] Our treatment is based on 
the quasiclassical theory of superconductivity. As in the 
vast majority of experiments the electron elastic mean 
free path £ in both superconducting and ferromagnetic 
layers remain much smaller than the coherence length 
^0 ~ vf/Tc, and ^fq ^ vp/h, respectively, where h is the 
exchange energy, it would be appropriate to stick to the 
so-called dirty limit and apply the Usadel equations^Si^. 
In addition, we will make use of the fact that at temper- 
atures sufficiently close to the critical one, T the 
superconducting order parameter A(T) is much smaller 
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Number 33 31 28 26 20 15 10 5 1 

dcuNi 2.5 3.4 6.6 10.8 23.9 39.6 60.1 79.5 89 

Tc (K) 6.73 6.7 5.25 2.403 2.686 2.855 2.975 3.005 3.349 

dH^2/dT (T/K) -0.5220 -0.4690 -0.4770 -0.2360 -0.216 -0.2210 - 0.2170 -0.2170 -0.2160 




FIG. 5: Temperature dependence of perpendicular (left panel) and parallel (right panel) critical magnetic fields for our FSF 
samples with the niobium layer thickness value dNb = 15.5 nm. With a good accuracy the thicknesses of the ferromagnetic 
layers are equal, i.e. di = d2 = dcuNi/2. Experimental data are shown by symbols. Solid lines represent fits to our theoretical 
predictions. 
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FIG. 6: (Color online) F1SF2 trilayer structure. 



than temperature, A(T) ^ T. In this limit, one can 
linearize the Usadel equations for the anomalous quasi- 
classical Green function F and obtain 



I? V-2i-A F 



2{\L0n\+ihsgnuJn)F + 2A{r) ^0, (1) 



where D = vp£/3 is the diffusion constant, A is the vec- 
tor potential, w„ = 7rT(2n -I- 1) is the Matsubara fre- 
quency and h is the exchange field in the ferromagnet. 
The exchange field h and the superconducting order pa- 
rameter A differ from zero, respectively, in the ferromag- 
netic and superconducting layers only. For simplicity, in 
what follows, we will assume that the magnetization in 
each ferromagnetic layer is spatially uniform within this 
layer, and that mutual orientation of the two magnetiza- 
tions is collinear, i.e. either parallel (P), or antiparallel 
(AP). 

Equation JT]) should be supplemented by appropriate 
boundary conditions matching the _F-functions at the 
S/F interfaces^. In the vicinity of the critical tempera- 
ture these boundary conditions read 



dF_ 
dx 



dF+ 
dx 



F. 



(2) 



where r is the interface resistance per unit area, a± and 
F± are the Drude conductivities and anomalous Green 
functions at both sides of the corresponding S/F inter- 
face, respectively. Below we will consider two orienta- 
tions of the external magnetic field: perpendicular and 
parallel to the S/F interfaces. 
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A. Perpendicular upper critical field 

Let us first consider an F1SF2 structure subject to an 
external magnetic field H perpendicular to the S/F in- 
terfaces. In order to account for a uniform magnetic field 
we can choose the gauge A = {0,0, Hy), where y, z are 
in-plane coordinates, and x-axis is perpendicular to the 
layers' plane. In this case we can use the following Ansatz 
for both the anomalous Green function and the order 
parameter— 



F(r,uJri) = g{y,z)F{x,u}n), 
A{r)= g{y,z) Aix), 

where the function g(y, z) obeys the equation 



(3) 
(4) 



d_ 

dz 



2i-Hy 

c 



4:Tt\H\ 



{m + l/2)g{y,z), m ^ 0, (5) 



and <I>o = TTc/e is the flux quantum. Then, Eq. ([ij takes 
the form 



D 



7^-^(1 + 2m) 



$0 



F 



2{\uJn\+ihsgi\uJn)F + 2A{x) =0. (6) 



Since the magnetic field enters only through the second 
term in Eq. (jG)), in order to obtain the critical field it 
suffices to restrict our analysis to the lowest Landau level, 
i.e. to set m = 0. 

In the ferromagnetic regions the superconducting order 
parameter A equals to zero identically, and Eq. ([6]) can 
be solved exactly. In the superconducting region it is 
necessary to obtain the solution of the equation 



dx"^ 



2t:\H\ 



Fs-2|w„|Fs + 2A(x) = (7) 



with the effective boundary conditions 
dFs 



dx 
dFs 



dx 



x=0 



x—d 



W^{u:)Fs{x^O), 



-W^{uj)Fs{x = d). 



(8) 
(9) 



where the parameters W^2 account for the influence of 
the top and the bottom ferromagnetic layers, respec- 
tively, on the superconducting layer. They read 



1 



1 



(Jo COthfci.2di.2 

ri,2 + 



1,2 



Cl,2 



Di, 



7tDi,2\H\ 



ihi,2 sgn w 



(10) 



(11) 



where (Ts^i.2 and I?s,i,2 denote the conductivities and 
the diffusion constants in the superconductor and in the 
two ferromagnetic layers with thicknesses d, di, and d2, 
respectively (see Fig. [S]). Exchange fields of the ferro- 
magnets are denoted by hi and /i2. The interface resis- 
tances per unit area ri and r2 refer to the S/Fi- and 
S/F2-interface resistances. 

Spatial dependence of the superconducting order pa- 
rameter is flxed by the self-consistency equation 



(12) 



\u)„\<u 



where A > is the BCS coupling constant, and uJc de- 
fines the high frequency cutoff, which is typically of the 
order of corresponding Debye frequency. Introducing the 
superconducting critical temperature Tco in the absence 
of both the external magnetic field and the ferromagnetic 
layers one can expel A and ujc from Eq. (|12p and get 



A{x)\n^^nT J] 



Aix) 



Fsix,UJn) 



(13) 



where T^q ^ uj ^{^-i j ■K)e-^ I ^ (7 = e"^' « 1.781). 

The critical magnetic field H;}^{T) is determined by 
a highest value of H at which the system of equations 
(H)-®, (fT3|) stfll has a non-trivial solution. In order 
to proceed further we will employ the fundamental so- 
lution method"^*. It is convenient to introduce the func- 
tion G{xi,X2,i^) which inside the superconductor (0 < 
xi,X2 < d) obeys the following equation 



dxf 



2 '^s 



G{xi,X2,uj) + d{xi - X2) ^ 0, (14) 

(15) 



2 firDslHl 



Ds V ^0 



together with the corresponding boundary conditions at 
the SF interfaces 



dG{xi,X2,u}) 



dxi 
dG{xx,X2,oj) 



2:1=0 



dxi 



WtG{0,x2,^), (16) 
^ -W^G{d,X2,u). (17) 



xi—d 



Equations (jl4p - (jl7p can be resolved analytically with the 
result 
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G{Xl,X2,Uj) = < r 



cosn/csxi + Wi(uj) ; 




coshks{x2 - d) - (cj) p 

ks 




ks 


sinh ksd + [W^{uj) + W^{uj)] cosh ksd 



cosh ks{xi 


^\ sinhfcs(a;i - 

-d)- W^{u;) 

ks 


d)] 




u 1 , w-L^ . sinhfcsa;2" 

coshfcs-T2 + W{-(lo) 

ks 






w^{uj)w^{ujy 

ks 


sinh ksd - 


f [W^^iuj) + W^{uj)] coshksd 



Xi < X2, 



X2 < Xi. 



(18) 



r 



The solution of Eqs. 0-® can now be expressed as a 
convolution of the function G{xi, X2,uj) and the spatiahy 
dependent order parameter A(x): 



which can be reduced by a partial summation over the 
Matsubara frequencies to a form: 



Fsix,io) 



G{x,x\uj)A{x')dx'. (19) 



Combining Eqs. ([T^ and we effectively reduce our 
problem to the integral equation for the superconducting 
order parameter which can be resolved in a straightfor- 
ward way. 

To proceed with numerical solution it is useful to dis- 
cretize this integral equation employing the Fourier trans- 
formation 



^(2;) = ^ \nfrn{x), 



fm{x) = COS 



m=0 

irmx 



< X < d. 



(20) 
(21) 



T ^ 



1 Ds\H\ 

2 2$nr 



i:Dsm\ 
ATd^ 



X A,„j(l +(5„ii,o) = X! 9mim2^m2, TOl = . . . OO, 

m2=0 



(23) 



Then, the integral equation can be transformed into an 
infinite system of linear equations, 



/mi ix)dx In 



TcQ 



\uj„\<(^c 

X / fmi{xi)G{xi,X2,UJn)fm2ix2)dXidX2 



where ip{x) is digamma function. The matrix elements 
9mi.m2 describe the strength of the proximity effect in the 
(22) superconducting layer. They are defined by the following 
expression, 



9mim2 



Ad^Ds 



[{ksdr + (^mi)2] [{ksdr + (^^2)2] 

[W^ + {-iy'^^+"'^Wi~]kssmhksd + W^W^[coshksd- + (-l)™i+™2] 

X ; ; — . (24) 



ks 



ks 



sinh ksd+ [W^ + W^] cosh ksd 



The critical magnetic field H^{T) coincides with the maximal value of the parameter H provided system of 
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equations (|23)) has a nontrivial solution. In our numeri- 
cal calculations an infinite system of equations (j23p was 
truncated to a finite one providing sufRcient accuracy of 
the solution. 



where xq is a free parameter which should be chosen in a 
way to maximize the critical magnetic field (or tempera- 
ture). Then, the Usadel equation ([T]) reduces to 



B. Parallel upper critical field 

Let us now turn to a configuration of the external 
magnetic field H oriented parallel to the S/F inter- 
faces. For a uniform magnetic field we choose the gauge 
A — {0,Hx,0) and employ the following Ansatz for the 
anomalous Green function and the order parameter-2^, 



D 



F 



-2{\ujn\+ihsgnuj„)F + 2A{x) ^ 0. (27) 



As before, we can effectively transform the above differ- 
ential equation into an integral one: 



F(r, uj) = exp ^— -^y F(x, lo), (25) 



A(r) = exp ( i ^^^^^ y ] A(a;), 



$0 



(26) 



Fs{x,uj) = — / G{x,x' ,ijj)l::^{x')dx' , (28) 





where now 



G(xi,X2,w) 



A2B2 qi{xi)qi(x2) 



AiBi 



q2{xi)q2 {X2) 



A1B2 - A2B1 



w 



+ < 



A1B2 - A2B1 w 
A2B1 qi(xi)q2{x2) , A1B2 q2{xi)qi{x2) 



A1B2 - A2B1 W A1B2 - A2B1 W 

A1B2 qi{xi)q2{x2) , A2B1 q2{xi)qi{x2) 



[ A1B2 - A2B1 



W 



A1B2 - A2B1 



w 



Xi < X2, 



Xi > X2, 



(29) 



r 



and 91.2(2^) denote two linearly independent solutions of 
the equation 



92 



dx^ \ <i>) 



f2nH\\ ,2 2\uj. 



q{x) = 0. (30) 



The parameters A\^2 and B\ 2 are defined by the follow- 
ing relations 



^1 =<7i(0)-Vt^J'9i(0), 
Bi = gi(rf)+W^I'Zi(rf), 
and 



^2 = <7^(0)-T4/llg2(0), (31) 
B2 = c{2{d) + Wlq2{d), (32) 



W = qi{x)q^{x) - q^{x)q2{x) 



(33) 



is the Wronskian for the two solutions q\(x) and q2{x). 
For simplicity we will assume that the orbital effects of 
the magnetic field in the ferromagnetic layers are small. 
Then, we can obtain parameters wf 2 from Eqs. ([TU]). 
([TT|) in which the external magnetic field should be set 
equal to zero H = 0. Combining Eq. ((29)) with the 
equations derived in Sec. 3A, we derive the upper critical 
field h}.2{T) for the parallel orientation of the magnetic 
field. 



On the basis of the analysis developed above the upper 
critical field Hc2 was evaluated numerically for the both 
cases of perpendicular (H^^iT)) and parallel (-ffiU^)) 
external magnetic fields. The corresponding results are 
presented in the next Section. 



IV. COMPARISON BETWEEN THE THEORY 
AND EXPERIMENT 



Let us perform a detailed comparison between our ex- 
perimental observations and theoretical results. By fit- 
ting our theoretical curves to the experimental data it be- 
comes possible to estimate the parameters of our F1SF2 
structures which could not be measured directly. Such 
fits for the perpendicular critical magnetic field measured 
in different F1SF2 samples are displayed in the left panel 
of Fig. [5] by solid lines along with the experimental data 
points denoted by symbols. The thickness of the S-layer 
remains the same for all samples in a series. The same 
set of parameters was employed to fit all the experimental 
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data presented in Fig. [5l These parameters are: 

Tco = 7.83K, ^s = 5.05nm, (34) 
hi^h2 = 99AK, (35) 
0-1 = era = 0.38f75, Di^D2 = i.2bDs, (36) 

Unfortunately, our fitting procedure does not allow to 
unambiguously determine the interface resistances ri^2- 
From Eq. ([T0|) we observe that under the condition 



< ri : 



< 



COth/ci 2^1 2 



0-1,2^1,2 



(37) 



the critical temperature remains almost insensitive to the 
values ri^2- Making use of the fact that the best fit of 
our data is obtained for fully transparent interfaces and 
employing the parameter values we arrive at an 

estimate 



< ri,2 ^ Cs/o-S- 



(38) 



Our theoretical results for the parallel critical magnetic 
field (solid lines) are presented in the right panel of Fig. 
[5] for the same values of fitting parameters 

In both cases of perpendicular and parallel magnetic 
field configurations we observe a reasonably good agree- 
ment between the theory and the experiment, in partic- 
ular, having in mind that single set of the fitting param- 
eters was used for all our samples. Further improvement 
of the quality of our fits can be achieved by a small ad- 
justment of these parameters individually for each sam- 
ple. We note that one can indeed expect such param- 
eters as exchange field /i, conductivity a and diffusion 
constant D to be slightly different for the ferromagnetic 
layers, since these parameters may depend, e.g.^ on the 
layer thickness and different growth conditions for the 
bottom and top Cu4iNi59 layers. There exists also a 
number of other physical reasons which might be respon- 
sible for a small mismatch between the theoretical curves 
and the corresponding experimental data in Fig. [S] One 
of such reasons is that the standard weak-coupling BCS 
theory of superconductivity employed here can describe 
the properties of niobium samples only approximatelySS: 
strong coupling corrections to this theory may easily 
reach ~ 30% in this case. This observation alone appears 
to be sufficient to account for the remaining discrepan- 
cies between the theory and experiment. In addition, 
the quasiclassical Usadel equation formalism employed 
here yields quantitatively correct results only provided 
the condition Ip ^ \J Dp/h is fulfilled in ferromagnetic 
metals. It is not completely clear if this condition is well 
satisfied in our samples, which are in the intermediate 
regime between strong and weak ferromagnets. Bear- 
ing all the above in mind we conclude that our theory 
sufficiently well describes the experimental data for the 
critical magnetic field in all our samples. 



To complete our analysis, in Fig. [7] we display both 
our theoretical results (solid lines) and the experimen- 
tal data (symbols) for the superconducting critical tem- 
perature Tc as a function of the total thickness of the 
ferromagnetic layers at different magnitudes of the ap- 
plied magnetic field. A good agreement between the the- 
ory and experiment is observed for all samples. At low 
magnetic field Tc (dcuNi) exhibits nonmonotonous depen- 
dence which becomes monotonous at high enough mag- 
netic fields. Similar dependencies are obtained for the 
critical field as a function of the total thickness of the 
ferromagnetic layers: a nonmonotonous dependence of 
^c2('^CuNi) at low temperatures turns to a monotonous 
one in the vicinity of Tcq. 

In conclusion, we have experimentally and theoreti- 
cally analyzed the upper critical field in F1SF2 trilayers 
in both perpendicular and parallel geometries. The series 
of samples with variable thicknesses of F-layers were pre- 
pared in a single run utilizing the wedge deposition tech- 
nique. Then, perpendicular H^{T) and parallel _ff].2(T) 
critical fields were measured in the same samples in the 
temperature range 0.4 — 8 K and for magnetic fields up 
to 9 Tesla. Temperature dependence of the critical mag- 
netic field is also derived from the quasiclassical Usadel 
equations for arbitrary values of the system parameters, 
such as superconducting transition temperature of the 
stand-alone superconducting layer Tco, superconducting 
coherence length ^5, exchange splitting energies hi 2 in 
the F-layers, conductivities of the ferromagnetic and su- 
perconducting layers cri,2,s, conduction diffusion coeffi- 
cients -Di,2.s in the F and S layers, and the interface 
resistances ri_2 for the S/Fi- and S/F2 interfaces, re- 
spectively. A good agreement between our experimen- 
tal data and theoretical predictions is demonstrated for 
all samples with variable thickness of the F-layers with 
the aid of a single set of fit parameters. The quasiclas- 
sical theory presented here allows to quantitatively de- 
scribe the experimental data for the upper critical field 
measured in F1SF2 hybrid samples alongside with Tc{d-p) 
and Tc((is) dependences and to collect an experimental 
database for determining a unique set of physical param- 
eters attributed to systems under consideration. 



Aknowledgements 

We would like to thank V.V. Ryazanov, 
M.Yu. Kupriyanov and A.I. Buzdin for useful dis- 
cussions as well as M. Schreck and S. Gsell for help 
with RBS measurements. The work is supported by 
A. von Humboldt Foundation within the frame of the 
Institutspartnerschaften Project and in part by the DGF 
grant No. GZ: HO 955/6-2. 



9 



\ ^\ \ ^ 


1 1 1 1 

H = 0T + 
i/ = 0.1Tx _ 
fl" = 0.2T « 
if = 0.3To - 
H = 0.4T ■ 


^ ' ■ 




^^^^^^^^^^^^^1^ — ET" 




J. + 


1 1 1 1 1 1 



15 20 25 




U = 0T + 
if = 0.6T X 
if = 1.2T « 
if = 1.8T D 
if = 2.4T ■ 



FIG. 7: Superconducting transition temperature versus the total thickness of ferromagnetic layers for both perpendicular (left 
panel) and parallel (right panel) magnetic field orientations. Different curves correspond to different values of the applied 
magnetic field (for clarity the curves corresponding to non-zero H are shifted upwards respectively by 2K, 4K, 6K and 8K). 
Experimental data points are denoted by symbols. Solid curves indicate the best fit to our theoretical predictions. 
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